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Abstract. The approach based on multimode system of g-deformed oscillators 
and the related picture of ideal gas of g-bosons enables to effectively describe 
the observed non-Bose type behaviour, in experiments on heavy-ion collisions, 
of the intercept (or the "strength") A of the two-particle correlation function of 
identical pions or kaons. In this paper we extend main results of that approach 
in the two aspects: first, we derive in explicit form the intercepts of n-particlc 
correlation functions in the case of q-Hose gas model and, second, provide their 
, explicit two-parameter (or qp-) generalization. 
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^} ' 1. Introduction 

o _ 

Quantum and g-deformed algebras are known to be very useful in diverse problems in 
many branches of mathematical physics and modern field theory 012] , as well as in 
O i' molecular/nuclear spectroscopy 3;- As well fruitful should be their direct application 

in the phenomenology of particle properties (see 0J0|B1IZ|j an d also [E] with references 
therein). Recently, it has been demonstrated [Hj that the use of multimode g-deformed 
oscillator algebras along with the related picture of ideal gas of g-bosons (q-Bose gas 
model) proves its efficiency in modelling the unusual properties of the intercept A of 
the two-particle correlation function, that is the measured value corresponding to zero 
relative momentum of two identical mesons, pions or kaons, produced and registered 
in relativistic heavy-ion collisions |1U) . where A exhibits sizable observed deviation 
from the naively expected purely Bose-Einstein type behaviour. The model predicts 
^DEL f° r a fixed value of q, the exact shape of dependence of the intercept A = A(K) 
on the pair mean momentum K and suggests asymptotic coincidence of A,r and Xk 
of pions and kaons. Put in another words, the intercept A, being connected directly 
and unambiguously with the deformation parameter q, tends in the limit of large pair 
mean momentum to a constant, lesser than unity, determined just by q and shared by 
pions and kaons. It is worth noting that confronting the predicted A,r behavior with 
the corresponding data from STAR/RHIC shows nice agreement ^2]> at least in the 
case of two-pion correlations. 

While two-particle correlations are known to carry information about the space- 
time structure and dynamics of the emitting source [TO] , in connection with some recent 
experiments it was pointed out 05] ^] that taking into consideration, in addition to 
the single particle spectra and two-particle correlations based analysis, the amount 
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of data concerning 3-particle correlations provides an important supplementary 
information on the properties of the emitting region, valuable for confronting 
theoretical models with concrete experimental data. Likewise, study of 4- and 5- 
particle correlations is also desirable [T5| . All that motivates the main goal of present 
contribution that is to derive the explicit formulas for the intercepts of higher order (n- 
particle, with n > 3) HBT correlations. Moreover, below we will obtain in explicit form 
the intercepts A*-™-* of n-particle correlations for the extended version of the developed 
approach when one uses the two-parameter gp-deformation of bosonic oscillators and, 
respectively, the model of gas of gp-bosons. 

From the very first, and up to more recent, applications of the g-algebras to 
phenomenology of hadrons there is growing evidence |SJ |5| that the phase-form of 
g-parameter is of great importance. Therefore, we hope that possessing the most 
general formulas for n-particle correlations and confronting them with the data from 
contemporary experiments will be helpful in clearing up the actual preference of 
choosing the form q = exp(i#) of deformation parameter. It is just this alternative 
for the choise of the deformation parameter q that implies very attractive physical 
interpretation of the g-parameter as the one that is directly linked to the mixing issue 
of elementary particles, either of bosons ^3E] or fermions |17l 1%]. 

The paper is organized as follows. Section 1 contains a sketch of necessary 
preliminaries concerning the two most popular types of g-deformed oscillators, as 
well as their two-parameter or gp-generalization. In section 2 we discuss basic points 
of the approach based on the g-Bose gas model, along with consideration of single 
particle g-distributions. The remaining two sections are devoted to the properties of 
two-particle and three-particle correlation functions, and to the main topic of present 
paper - the results on the multi-particle (n-th order) correlations, for the algebras of 
both the g-deformed and the gp-deformed versions of generalized oscillators. Details 
of derivation of basic formulas are relegated to the Appendix. 

2. g-Deformed and gp-deformed oscillators 

We begin with a necessary setup concerning two types of g-deformed oscillators, and 
also their two-parameter generalization. 

g-Oscillators of AC type. The g-oscillators of the Arik-Coon (or AC-) type are 
defined bv the relations 18! FT9j 

a^a] - q 5,j a] j a i = Sij [ai,aj] = [a\, at] = 

Wi, = Wi,a]} = Sijd] Wi,Nj] = . (1) 

were — 1 < q < 1. Note that this is the system of independent g-oscillators as clearly 
seen at i =/= j. 

From the vacuum state given by a^jO, 0, . . .) = for all i, the basis state vectors 

\m,...,m,...) = 1 (a\r (4r • • • (4r • • • 10, o, . . .> (2) 

V [nij ! L"2j ! • • ■ \nij\- ■ ■ 
are constructed as usual, so that 

a| I ...,m,...) = VT^i + lJl • • • ,n t + l, ...) (n\ . . . ,ni, . . .) = \f\jH\\ ■ ■ .,n<-l, . . .) 
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Mi\n l7 ...,n,i,...)= rii\ni, ...,m,...) (3) 

Here the notation [• ■ -J for so-called basic numbers and the corresponding extension 
of factorial, namely 

W = y=7 LrJ! = LiJL2j---lr-iJlrJ Loj ! = LiJ ' = i (4) 

are used. The g-bracket [A\ for an operator A is understood as formal series. At 
q — ► 1, from [^J and [A\ one recovers r and A, thus going back to the formulas for the 
standard bosonic oscillator. For the deformation parameter q such that — 1 < q < 1, 
the operators a\ , at are conjugates of each other. 

For g 7^ 1, the bilinear a\ai depends nonlinear ly on the number operator J\fi : 

a\a t = LMJ (5) 

so that at q = 1 the familiar equality a\ai — Mi is recovered. 

^-Oscillators of BM type. The g-oscillators of Biedenharn-Macfarlane (BM) type 
are defined by the relations |2UM19) : 

[h, b 3 ] = [bib]} = [N u b 3 ] = -b l3 b, [Ni, 6t] = fyftt [N u Nj] = 

bib] - q S * b]bi = Syq-"' hM\ - q^ b]b t = 6^ . (6) 

In this case the extended Fock space of basis state vectors is constructed in the way 
similar to the above case, with the only modification that now we use, instead of basic 
numbers, the g-bracket and g-numbers, namely 

b\bi = [N& [r] q ee . (7) 

Formulas similar to (2)-(5) are valid for the operators bi, b] if, instead of (4), we now 

use the definition (7) for g-bracket. Clearly, the equality b\bi = Ni holds only in the 
"no-deformation" limit of q = 1. For consistency of the conjugation, we put 

q = exp[i6) < 9 < 7T . (8) 

gp-Oscillators. Besides the g-bosons of AC-type and BM-type, in what follows 
we will also consider the two-parameter (or qp-) generalization of deformed oscillators 
given by the relations [21] 

[JV(«0 , A] = —A [N { ^ , A 1 "] = A* 

AA* - qA^A = p N AA f - pA^A = q N (9) 

from which 

x x 

At A = lN^j qp with \X\ qp ee q —- *— . (10) 

In this definition we have shown only one mode, although similarly to (1) and (6), 
in what follows we will deal with the system of independent (that is, mutually 
commuting) copies/modes of the gp-deformed oscillator. Note that X in (10) can 
be either a number or an operator. Clearly, putting p = 1 immediately leads us to the 
AC-case while putting p = g -1 reduces to the BM-type of g-bosons. 
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3. Statistical (/-distributions 

For the dynamical multiparticle (say, multi-pion or multi-kaon) system, we consider 
the model of ideal gas of q-bosons (IQBG) by taking the free, or non-interacting, 
Hamiltonian in the form pfl 1321 EH] 

H = J2"iMi (11) 



where uji = y/m 2 +k 2 , A/i is the number operator given in (5) or (7) or (10), and 
the subscript V labels different modes. Let us note that among a variety of possible 
choices of Hamiltonians, the choice (11) is the unique truly non-interacting one, which 
possesses an additive spectrum, see |22I QH] ■ Clearly, it is assumed that 3- momenta of 
particles take their values from a discrete set (i.e. the system is contained in a large 
finite box of volume ~ L 3 ). 

To obtain basic statistical properties, one evaluates thermal averages 

where (3 — 1/T and the Boltzmann constant is set equal to 1. Calculating, say, in the 
case of AC- type g-bosons the thermal average (q^* ) , with A/j from (5) , with respect 
to the chosen Hamiltonian (11) we obtain 

„r e^' - 1 

<^'> - < 12 > 

and the distribution function (recall that —1 < q < 1) is found as |22l I23| : 

Usual Bose-Einstein distribution corresponds to the no-deformation limit of q — > 1. In 
the particular cases q = — 1 or q = the distribution function (13) yields respectively 
Fermi-Dirac or classical Boltzmann ones. Note that this coincidence is rather a formal 
one: the defining relations (1) at q = — 1 or q = differ from those for the system 
of fermions or the non-quantal (classical) system. The formal coincidence of equation 
(13) at q = —1 with the Fermi-Dirac distribution can be interpreted in terms of 
impenetrability (or hard-core) property of such bosons. The difference with the system 
of genuine fermions lies in commuting (versus truly fermionic anticommuting) of non- 
coinciding modes at q = — 1, see (1). 

Now consider BM-type of q-bosons. The Hamiltonian is chosen again as that of 
IQBG, but now with the number operator given in (7), i.e., 

H = J2"iNi- (14) 

i 

Calculation of (q ±N >) yields {q ±N *} = (e^ - l)(e? u < - ? ±1 )" 1 . Then, from the 
formula {b\bi) = (eP u)i — q)^ 1 (q~ Ni ) the expression for the (/-deformed distribution 
function (note that q + q^ 1 = [2] q = 2 cos 8) does follow (see also |22l 123) ): 

{bih) = e 2 ^-2cos0e^ + l ' (15) 
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Although the deformation parameter q is taken as complex one according to (8), the 
explicit expression (15) for the q-distribution function shows that it is real. 

It is easily seen that the shape of the function /(k) = (b^b)(k) in (15) is such 
that the q-deformed distribution function with q ^ 1 is intermediate relative to the 
other two curves, the standard Bose-Einstein distribution function and the classical 
Boltzmann one (the same is also evident for the above q-distribution function (13) 
of the AC-type q-bosons). That is, the deviation of the q-distribution (15) from the 
quantum Bose-Einstein distribution goes, when q goes away from the no-deformation 
limit q = 1, in the "right direction", towards the classical Boltzmann one. 



4. Two- and three-particle correlations of q-bosons 



Although the formulas for two-particle correlation functions have been obtained earlier 
[HI, we recall them here for the sake of more complete exposition. In the remaining 
part of this section some new results will be presented. So, we consider two-particle 
correlations first with the AC- type of q-bosons. Starting with the identity 

a\a]akai — q~ 5ih ~ 5il a]akaia\ = [a] , a]]a fc a; + at[a|, ak] q -s ik ai + q~ 5zk a^a k [a\, ai] q -s u 

where [X, Y] K = XY — kYX , by taking thermal averages of its both sides we find 

gfr* - ■ 
q^.-u ik -5u e ^i — q 

For coinciding modes this leads to the formula 



{a\a]a k ai) = ^JLg,, ^ z {{a]a{) (a\a k ) + q s ^ (a]a k ) (a\ai)) 



Mm) - (e ^_ g) | e L- g2) ■ (16) 

From the last relation and the g-distribution (13) the ratio under question (called 
intercept) does result: 



(fljjojflj) 1= 1 | (l + g)(efr* -g) _ e^ - 1 
(afa,) 2 " e^-q 2 - q 2 



Note that in the non-deformed limit q — * 1 the value Abe = 1, proper for Bose-Einstein 
statistics, is correctly reproduced from equation (17). This obviously corresponds to 
the Bose-Einstein distribution contained in (13) at q — ► 1. The quantity (intercept) A 
is important since it can be directly confronted with empirical data. In this respect, 
let us note that there exists a direct asymptotic relation A = q, which corresponds to 
the limit of large momentum or low temperature (in that case f3oj — > oc). 

Now we go over to the Biedenharn-Macfarlane q-oscillators (6) and find the 
formula for the monomode two-particle correlations, i.e. for identical particles with 
coinciding momenta. From the relation 

(bjbtbA) - q 2 {b\bAb\) = -{b\b iq ^){\ + q 2 ) 
valid for the monomode case at hand, we deduce 
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Evaluation of the thermal average in the r.h.s. yields (b\biq Ni ) = q/(e^ UJi — q 2 ). Using 
this we find the expression for two-particle distribution, namely 

;,t,+, , i 2cos# 

= e ^-2cos(2^ + l • < 18 > 

Then, the desired formula for the intercept of two-particle correlations of the BM-type 
g-bosons, with the notation ti = cosh(/3cji) — 1, reads 

x 1 , (b\b\bA) 2cos^ + l-cosg) 2 
{(blh)? *? + 2(1-008^ 

and again is a real function. 

Three-particle correlations of the q-bosons of AC-type 

Derivation of three-particle correlation functions proceeds analogously to the 2- 
particle case. Considering the g-deformed oscillators of AC-type we start with the 
easily verifiable identity 

a^alaia m a s al = a { j a\aia m [a s ,a { i \ q s zs + 
+q s " {alal ai ([a m ,al] q s tm )a s + q Sim {a]a\{[ai, a\\ qSil )a m a s + q 5u a]a\a\aia m a^ j 



and take thermal averages of both its sides. This leads to the equality 

e iJ^i _ q 



i t t t \ e> z - q /, t t wt\ 

{a l a J a k aia m a s ) = _ qSis+Stm+Su \{ay k aia m ){a\a s ) 

+q s " (a]a\aia s )(a\a m ) + q 5is+5im (a)a\a m a s ){a\ai)} 



which in view of (a\(ij) — 5ij(a\ai) — S i j/(e l3LJi — q), cf. (13), in the monomode 
i=j = k = l = m = s case yields: 

/ t t t v (l + g)(l + g + g 2 ) 

(alalala^at) = (e ^, _ g)( ^, _ g2)(e ^ _ g3) ■ ( 20 ) 

From the latter relation, dividing it by (a\ai) 3 , we derive the desired expression for 
the intercept (or strength) A^ 3 ' of 3-particlc correlation function (we drop the V): 



A (3)_ (a t3 a 3 ) , (l + q)(l + q + q 2 )(e^- q r 1 



In a similar manner, it is possible to derive the (intercept of) 3-particle correlation 
function for the system of BM-type g-bosons. However, intead of doing this, in the 
next section we will derive the most general results for both 3- and n-particle, n > 3, 
correlation functions in the two-parameter (i.e. gp-deformed) extension of bosons, 
from which the desired formulae for the BM-type of g-bosons will follow as particular 
cases. 
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As extension of equations (16), (20), it is not difficult to derive, using method 
of induction, the following general result for the n-particle monomode distribution 
functions of AC-type g-Bose gas: 

((4) n ( a i) n ) - n^^-gr) H = = i + q + f + ... + r - . ( 22) 

From this expression the desired formula for the intercepts = ^^pr — 1 of 
n-particle correlations of AC-type g-bosons immediately follows (with V dropped): 



n"= 2 (e^-9 r ) 



= ^ ~^ ■ (23) 



In the asymptotics of /3w — > oo (i.e., for very large momenta or, at fixed momentum, 
for very low temperature) the result depends only on the deformation parameter: 



A H asympt = -1 + [„j ! = -1 + JJ 

k=l V=0 



k 

V 



= (1 + g)(l + q + q 2 ) ■ ■ ■ (1 + q + . . . + g™" 1 ) - 1. (24) 

This remarkable fact can serve as the test one when confronting the developed 
approach with the numerical data for pions and kaons extracted from the experiments 
on relativistic heavy ion collisions. 

Now we come to the base point. 

Extension to qp-bosons. 

The above results admit direct extension to the case of the two-parameter 
deformed (or qp-) oscillators and thus to the gp-Bose gas model. For this, we use 
in analogy with (11) and (14) the Hamiltonian 

H = J2^N^. (25) 

i 

With (25), the expression for general n-particle distribution functions is obtained (see 
Appendix for its derivation) as 

q rn _ p m 

\m\ qp = _^ \m\ qp \ = [l] 9 p|2]gp • • • [m - l]gp[m] 5P . 

In the particular cases n = 1 and n = 2 (note that [2] ?p = P + q) this obviously yields 
the formulas 

i e fai _ i) 
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i U2 / A^_ (p + g) (e*"<-l) 



( e /3^ s _ qt)(eP<n -pq)( e P«>i - p 2 ) 



(remark that the latter two formulas were also found in |25|1. 

From (26), after dividing it by {A\Ai) n , the most general result for the n-th order 
<7p-deformed extension of the intercept \( n \ omitting the V, follows as 

x (n) _ ^ 2 ' - 1 = M ' - P> - 1 (27) 

q,p ~ (a^a) u l lqp -( e ^-iy^u r k=o( el3uJ -i n ~ k p k) 



which constitutes our main result. This provides generalization not only to the case 
of 7i-th order correlations but also to the two-parameter (gp-)deformation. 

Let us give the asymptotical form of intercepts in this most general case, X q n p'. 

n / k \ 

A („), asympt = _j_ + ^ = -1 + JJ ( £ jV" ) . (28) 

fc=l V=0 ' 

As we see, for each case of deformed bosons (the AC-type, the BM-type, and their 
gp-generalization) the asymptotics of the n-th order intercept takes the form of the 
corresponding generalization of the usual n-factorial (the latter yields pure Bose- 
Einstein n-particle correlation intercept). 

Finally, let us specialize the obtained formulae to the case of g-bosons of BM type 
for n = 3, that is 

A (3) = , |A (e 2 ^-2e^cosfl + l) 2 

(eP» - l) 2 (e 2 ^ - 2e<3" cos(3#) + 1) [ ' 

A (3), asympt = _i + [ 2 ] g [3] 9 = — 1 + 2 cos 6 (2 cos - 1) (2 COS 9 + 1).(30) 

In conclusion we note that it would be of great interest and importance to make a 
detailed comparative analysis of the obtained results with the existing data for 3- 
particle correlations of pions and kaons produced and registered in the experiments 
on relativistic heavy ion collisions, with the goal to draw some implications concerning 
immediate physical meaning and admissible values of the deformation parameters p, q. 
Details of such analysis will be presented elsewhere. 

Appendix 

Here we derive the general formula, see (26), for the (monomode) n-particle pq-bosonic 
distribution functions: 

<„...„.) . . (A „ 

For convenience, in (A.l) and below, we drop the mode-labelling subscript "i", and 
use a\a,N instead of A\ A, N^ qp ^ respectively. The proof proceeds in few steps. 
First let us derive the recursion relation 

(a+"a«) = (a^-V-y) , „ Wgp N T . (A.2) 
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For this, we use pg-deformed commutation relations and evaluate the thermal averages: 

= (a^-Wa"- 1 )^ - (a^ n - 1 a n - 1 p N }^— 
x 1 q x F 'q-p™" 1 

= (a*" Wa»" a )l - - f 4i + -k) (« t - 1 «"-V Ar ) = • • • 

g ^p™ 1 qp™ 1 J 

= (a*"- W)- - I f J_ + _J_ + ... + (at"- 1 ^"^) 
= < a t»-W) — - ( a t»-V-i^)_N«L 

= (a tn a n >^ - (« t "" 1 a""V JV >^fe • (A.3) 

From this the equation (A. 2) readily follows. After k-th iteration of this procedure we 
find 

{a tn-k a n-k k N) = (a tn-(fc+l) a n-(fc+l) (fc+l)JVv [ " 7 ^ ^ ■ (A.4) 

> -r- / \ -r / /g/Sa) qn—kpk\pn — (2k+l) v y 

Indeed, 

(at"- fe a"- fe p feJV ) = (a tn - (fe+1) aat a "- fe -V feAr )i - ( a tn-(k+i) p " a »-(fc+i)pfcJV)i = ... 

P 0u 1/1 



+ - ^ / a tn-(fe+l) a n-(fe+iyfe+l)W\ 

n — h — 1 J 



= (at-'a-V")^ - (^- (fc+1) ^- (fc+1) ^ (fc+1)JV ) ^%T) (A.5) 

that is equivalent to the formula (A.4). Applying this formula step by step n times 
yields the relation 

(at"a") = = ^ = ( P nN ) . (A.6) 

nVo(^-p r <i n - r )m=oP n - (2k+l) 

From the latter, with the account of 

(P nN ) = fiV-^ - 1 (A.7) 

e p k=0 

wc finally arrive at the desired formula (A.l) for the higher order (n-particle) 
monomode distribution functions of the model of gp-Bose gas. 
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